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This paper is based on results derived during a stability study of the 
Saturn V rocket jor ivhich it ivas necessary to validate the use of Nyquist's 
encirclement-counting technique in distributed systems. An outli?ie of the 
paper is as follows: Certain results concerning the finiteness of the number 
of zeros of polynomials in s and e" are shown in Theorem 1 and its cor- 
ollaries. Theorem 2 is a generalization of Michailov's Criterion. Simplifying 
assumptions, usually valid in practice, yield a simplified test to determine 
if "encirclement-counting" is a valid stability test [equation (20)\. The 
results are reformulated for an open-loop analysis. Various aspects of the 
theory are shown by three examples based on an electrical equivalent of a 
simple single-engine, liquid-fuel rocket. 

I. INTRODUCTION 

Liquid-fueled rockets can exhibit a peculiar type of instability due to 
self-sustained longitudinal oscillations. Since the rocket then stretches 
and shrinks longitudinally, it behaves like a pogo-stick, which has 
resulted in the nickname POGO for this type of instability. 

To see how this phenomenon arises, consider the simple diagram 
shown in Fig. 1. The chain of events which can cause POGO is initiated 
by a random variation in thrust of the engine. This thrust variation 
causes the rocket structure to oscillate in its natural modes. The pressure 
in the fuel tank thus varies. This pressure variation is propagated down 
the fuel feed line, resulting in a variation of fuel flow into the engine. 
Since the thrust of the engine is proportional to the rate of fuel entering, 
the loop is completed, and instability results if this resulting thrust 
variation aids the original random variation which initiated the chain 
of events. 



* This work was done when the author was with Bellcomm, Inc., and therefore 
was performed under NASA contract NSW417. 
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ROCKET 



FUEL TANK 



FUEL FEED LINE 



/ V- ENGINE 

Fig. 1 — Simple model of a liquid-fuel rocket. 



In modeling this physical situation, it is customary to perform a modal 
analysis of the rocket structure, and retain only the most significant 
modes. This results in the distributed structure being replaced by a 
lumped approximation. Although the same technique could be used 
to lump the feed-line, there are many reasons for desiring to keep this 
element as a distributed parameter. There is no great difficulty in doing 
so, since the fluid equations which govern the feed-line are of the same 
form as electrical transmission lines. 

An equivalent circuit of Fig. 1 would then be Fig. 2, where Vi(s), V 2 (s) 
are the Laplace transforms of the pressure variations at the top and 
bottom of the feed-line, respectively, and Ii(s), / 2 (s) represent the 
transform of flow variations. EJs) is then the random pressure variation 
at the top of the line due to the assumed random thrust variation above. 
Y(s) represents the hydro-mechanical impedance of the feed-line output, 
and G(s) includes the structural feedback. To see the type of equations 
which will be of concern for a stability analysis, assume that losses can 
be neglected in the feed-line. 
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Fig. 2 — Electrical equivalent of Fig. 1. 
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The system equations can be written in matrix form as: 
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(1) 



where Z c (a positive real number) is the characteristic impedance of the 
line and k (a positive real number) is the ratio of line length to wave 
speed in the line. The determinant (A(s)) of the matrix appearing in (1) 
is of interest in a stability analysis. It can easily be computed to be 



A(s) = cosh ks - fr'(.s) 4- Z r Y(s) sinh ks. 



(2) 



In general, neither Y(s) nor G(s) need be rational (especially in the case 
of multiple engines) but for simplicity, suppose that 



G(s) = 



Ng(s) 

D G (s) 



Z c Y(s) = 



N Y (S) 

Dy( S ) 



(3) 



where N r , , D a , N Y , D Y are polynomials in the complex variable s. 
Then (2) can be written as 

A(s)[D G (s)D Y (s)} = [D a (s)Dy(s)] cosh ks - N G (s)D y (s) 

+ D G (s)N Y (s) sinh ks. (4) 

In the following sections of this paper we will be concerned with 
polynomials in the two complex variable s and e '. We now show that (4) 
can be considered as such. Thus multiply (4) by 2e _2t " to give 

F(s) = 2e- k ' A(s)[D G (s) D y (s)] = D G (s)[D Y (s) 4- N Y (s)] 

- 2N G (s) D Y (s)e- k! + D G (s)[D Y (s) - N Y (s\]e- 2k ' (5) 

F(s) = Ei2,(s>-*' (6) 



or 



where each R ( (s) is a polynomial in s. 

This paper will also be concerned with open-loop and closed-loop 
expressions, and will assume that the quantities of interest will be of 
the form of (0). To show that this is true in our present example, assume 
that G(s) "closes the loop" and solve for the "line transfer function" 
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7 a (s)/V,(s). The result is 

7,(«) 1 



Vi(s) cosh ks + Z c Y(s) sinh ks 
Hence the "loop gain," G n (s), is simply G(s) times (7), or 

m 



(7) 



G (s) = 

(8) 



D Y (s) D (s) cosh fcs 4- D (s)N Y (s) sinh A;s 

The denominator of (8) is seen to be of the form of A(s), and hence can 
be made to look like (6). 

Finally, we wish to remove the simplifying assumption that Y(s) 
and G(s) are rational. This is desired since an actual engine is not only 
fed liquid fuel, but also liquid oxidizer. Thus even a single-engine rocket 
has two feed-lines. (The Saturn V has five engines, for a total of ten 
feed-lines.) The structure of a model for rockets of this complexity 
would be that of Fig. 2, repeated once for each feed-line, with suitable 
interconnections through lumped (i.e., rational) transfer functions. 
Thus the G(s) and Y(s) of Fig. 2 will be the ratio of sums of powers of s 
and e". As such they can be combined to yield forms such as (6). 

A final constraint on the stability analysis is that it is required that 
the analysis be of the conventional open-loop type using Nyquist's 
Criterion. For lumped systems this presents little difficulty since one can 
always make open-loop measurements at as high a frequency as neces- 
sary to guarantee that all singularities of the transfer function are 
included. Furthermore, for rational functions, no difficulty is encountered 
in closing the contour in the right-half s-plane. With distributed systems, 
however, it is possible that the gain becomes periodic for large magnitude 
of s, and care must be exercised in determining closed-loop stability 
via open-loop gain plots. What is desired, therefore, is a set of conditions 
under which the conventional "encirclement counting" technique for 
lumped systems, remains valid for distributed systems of the type 
described. 

There are two important techniques for determining whether a 
polynomial in the two complex variables s and e* has any zeros for 
Re [s] ^ 0. These are the Pontryagin Criterion 1 and the Michailov 
Criterion. 2 It is of interest to see if these criteria can be applied to the 
ratio of such polynomials in order to determine stability of closed-loop 
gain. In a previous note 3 it was shown that this is not feasible for criteria 
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of the Pontryagin type. In this paper we are able to develop a stability 
criterion of the Nyquist type from Michailov's Criterion for a large 
class of distributed-parameter systems, in particular, for a large class of 
transmission line systems with feedback. 

An outline of the paper is as follows: Certain results concerning the 
finiteness of the number of zeros of polynomials in s and e are stated 
as Theorem 1 and its corollaries. Theorem 2 is the desired generalization 
of Michailov's Criterion. Simplifying assumptions, usually valid in 
practice, yield a simplified test to determine if "encirclement-counting" 
is a valid stability test [equation (20)]. The results are reformulated 
for an open-loop analysis. Three examples show various aspects of the 
theoretical analysis. The Appendix includes a statement of Pontryagin's 
Criterion suitable for use in the present paper. Also included are the 
proofs of the two theorems and a derivation of some conditions under 
which Michailov's Criterion can be simplified. 

ii. michailov's criterion 

As a starting point we consider an equation of the form 

G(z)= EEfliA"" = (9) 

1=1 J=0 

where a,, are complex and w, are real. If any of the co, were negative, 
we could multiply G(z) by e'" fcU , where w t is the most negative of the cos. 
This would not change the zeros of G(z), so we assume ^ to, < w 2 < 
• • • < co m . Dividing by e Um ' and letting a m _ 1 + 1 ,, = a,, , we transform (9) 
into 

m = E i; a t ^e-'" (10) 

1=1 1=0 

where r, = to m — w m _ 1+1 > for i = 2, 3, • • • , m and r x = 0. [To relate 
this to the Pontryagin Criterion, note that if the co, are rational (which 
can always be assumed in a practical situation) then a suitable scaling 
of the z variable will make G(z) of (9) into a polynomial like H(z) of the 
Pontryagin Criterion.]* 

Before continuing, it should be noted that a proof of the Michailov 
Criterion for exponential polynomials has been presented in the litera- 
ture. 2 However, this proof assumed that 

I a,. | > E|a,„ |. (11) 



* See Appendix and Ref . 1 . 
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As will be seen by the example to be considered later, in transmission 
line systems, (11) is almost never satisfied. Hence, a proof is desired 
which is free of this assumption. However, we do use the assumption 
that Oi„ 5^ 0. That this is no loss of generality can be seen by the follow- 
ing considerations. It is clear that | r m | Si | r t | , i = 2, • • • , m. Then 
multiplying (10) by e mZ (which does not change the zeros) puts (10) 
into the Pontryagin form. If a ln is zero, the principal term is missing 
and we are finished with the stability study.* To aid in the subsequent 
development, let us rewrite (10) as 

F(z) = Z^.n (12) 

1=0 

where Qi(e~ z ) = ^?=i a 1 -,e~ r ' a , or as 

F(z) = z k b\{z) + E *'Q,(0 (13) 

where k is an integer between zero and n and 

F k (z) = £y-*Q f -(e"). (14) 

i = k 

In the Appendix we prove the following theorem: 

Theorem 1: If there exists a non-negative integer k ^ n such that F k (z) 
of (14) has at most a finite number of zeros on Re (z) 2: 0, then F(z) of (13) 
has at most a finite number of zeros on Re (z) > 0. 

The following corollaries are of interest: 

Corollary 1: If there exists a non-negative integer k ^ n such that F k (z) 
of (14) has no zeros on Re (z) 2: 0, then F(z) of (12) has at most a finite 
number of zeros on Re (z) > 0. 

Corollary 2: If Q„(e~') of (12) has no zeros on Re (z) ^ then F(z) 
of (12) has at most a finite number of zeros on Re (z) > 0. 

[Corollary 2 follows since Q n (e~ z ) — F n (z).] 
Let 

F(z) = 1 + +(z)* k F k (z) (15) 

where 

*M -£«"-"%£?• (16) 

,=o tk\t) 

In the Appendix we derive theorem 2: 



* See Appendix and Ref. 1. 
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Theorem 2: The number of zeros of F(z) with positive real part is 

+ ± arg (1 + Wu)) ~ arg (1 + *(-*y)) (17) 

assuming that F(z) has no purely imaginary zeros and where b. u (F(z)) 
is the net change in arg F(z) along the imaginary axis from — iy to +iy, 
and A c (F k (z)) is the net change in arg (F k (z)) along contuor C, any contour 
outside the semicircle of radius R of Theorem 1. 

III. MICHAILOV'S CRITERION! SPECIAL CASE 

Theorem 2 is the desired statement of Michailov's Criterion. To 
obtain tighter results, let us now assume that Q„(e~') has no zeros on 
Re (z) ^ (i.e., Corollary 2). Further, let the r, be rational and the 
a,,- be real. By virtue of rational r, , Q k (e~') is periodic in y, for z = 
x + iy. Let this period be P. By virtue of real a,, , replacing z by its 
conjugate results in Q k (e~') being replaced by its conjugate. Hence we 
need only consider the semi-infinite strip defined by x > and P ^ y ^ 0. 

Michailov's Criterion can be simplified if it can be shown that Q n (e~') 
does not wind around the origin as z varies over a suitable C. We now 
consider this possibility. 

Q n (0 = E a,„e~ r "(cos r,y - z'sinr,?/) 
i-i 

m n» 

= a ln + JJ a in e" riz cos r t y - i J] a in e~" z sin r,y 

,=2 i-2 

= Re[Q n ] - ilmLQJ. (18) 

If either Re [Q„] or Im [Q„] does not vanish along C, then Q n cannot wind 
around the origin. In the Appendix we derive sufficient conditions for 
this. 

We now assume that A r [Q n (e~ : )] = and write the Michailov Criterion 

as 

N = I - - A_„, /2 (F(*)) + ~ arg (1 + f(iy)) (19) 

A IT IT 

where A_„, /2 is that part of the imaginary axis from to iy. Thus N = if 
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and only if 

A. w/2 (F(z)) = *f + arg (1 + *(*/)). (20) 

We remark that arg (1 + $(iy)) can be made close to zero for y suf- 
ficiently large. 

IV. MICHAILOV'S CRITERION APPLIED TO OPEN-LOOP ANALYSIS 

The Michailov Criterion, as well as its predecessor the Pontryagin 
Criterion, settle the problem of finding rhp zeros of polynomials in z 
and e 1 '. In many engineering applications, however, this polynomial is 
not directly available, but a related ratio of such polynomials can be 
found. In the study of feedback systems, for example, an open-loop 
gain can be measured and it is desired to find the poles of the closed-loop 
gain. These latter poles are the zeros of the polynomial which results 
from adding the two polynomials whose ratio is the open-loop gain. 
Stability has been determined for nondistributed systems by counting 
encirclements of the open-loop gain along the imaginary axis. What we 
propose to do next is to provide a similar criterion for the distributed 
parameter problem. Thus let F{z) = D(z) + N(z). Then 



A r (D(z) + N(z)) = A r (Z)( 2 ))(l + ^ 

= A r (Z>(*)) + A r (l+fj|)- (21) 

Let the contour r be composed of a portion of the imaginary axis w 
and another (possibly semicircular) contour C, such that T encloses 
all zeros of D{z) + N{z). Then 

A r (D(z) + N(z)) = A r (D(z)) + A„(l + |g) + A c (l + |®)- (22) 



It is the term 



mi -m 



D{z), 
which is usually available for determining stability. We ask, 

"When does A w {\ + ^|j) = A r (D(z) + N(z))V. 
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The answer is that this happens exactly when 

= A r (D(*)) + Ac(l+|^j)- (23) 

To develop a more practical criterion let us rewrite this expression 
using n F to be the highest power of z in F(z), and N F to be the number 
of zeros of F(z) inside r. Then 

= A r (D(*)) + Ac(D(*) + N(z)) - A c (D(z)) (24) 

= 2tN d -\- n D+N ir — n D ir (25) 

where we have neglected those terms which become small for large z. 
If we limit further consideration to systems which are open-loop stable 
(i.e., N D = 0) then (25) requires that n D+N = n D . In most practical 
situations, the open-loop gain is bounded at infinity, that is to say 
n N 1= n D . Hence n D+N ^ n D . Since n D+N < n D requires 

llm DM = _1 > 

we can conclude that counting encirclements of the open-loop gain is a 
valid method for determining stability [i.e., (25) is satisfied] for systems 
which are open-loop stable (N D = 0) and whose gain is bounded at 
infinity (n D+N ^ n D ) but does not approach —1 for large frequencies 
(n D+N < n D ). This includes the case, usually found in practice, that the 
open-loop gain approaches zero for large frequencies. 

v. EXAMPLES 

All examples refer to Fig. 2 and are chosen to illustrate various aspects 
of the analysis. At various points in the examples, the following assump- 
tions concerning G(s) and Y(s) are referred to: 
Al. G(s) and Y(s) are each the ratio of two polynomials having real 

coefficients with no singularities on Re (s) > 0. 
A2. lim.-Hjo G(s) = ki where /c t is real and | k v \ ^ 1. 
A3, lim,-,,,, Y(s) = lim.-.a, sC where C is non-negative real. 

Assumption Al requires G and Y to be stable transfer functions. 
Assumption A2 insists that the feedback gain at infinity be less than 
unity. Assumption A3 is physically appealing. 

Some unusual properties of the natural frequencies of this system 
have been described elsewhere. 4,5 

Exam-pie 1: In this example, Assumptions Al, A2, and A3 are invoked, 
and G(s) and Y(s) are given by (3). We wish to show that assumption 
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(11), used in previous derivations of Michailov's Criterion, is not met 
and that it is valid to count encirclements of the Nyquist plot to de- 
termine stability. By Al, D (s) and D Y (s) have no right-half plane 
zeros. Thus A(s) has right-half plane zeros exactly when A(s)D (s)D Y (s) 
has right-half plane zeros. 

From Assumption A2 we conclude that deg D (s) ^ deg N a (s). 
From A3 we conclude that deg N Y (s) > deg D Y (s). From this, and (5) 
and (6), we see that the principal term is present and that the assump- 
tion (11) used in previous proofs of Michailov's Criterion is not met. 
In fact, one can readily convince oneself that this will be the case when- 
ever lossless transmission lines are involved, since all exponential terms 
will involve hyperbolic functions. 

Using the notation of (6), the open-loop gain (8) can be expressed as 

G (s) = o rffirL-^ < 26 ) 

R {s) + R 2 (s)e 

whose norm becomes small for large, right-half plane values of s. 
Hence it is valid to count encirclements of the open-loop gain about 
the point +1. 

Example 2: Here we show the necessity of the open-loop stability 
requirement. Suppose Assumptions Al and A2 are invoked and further 
assume that D a (s) = 1, N G (s) = k t , D Y (s) = 1, N Y (s) = sC x + g. 
This corresponds to terminating the line in a capacitance Ci and shunt 
conductance g ^ 0. The open-loop gain for s = jw becomes 

PC \ f^i C27) 

°^ ' ~~ cos kw — wCi sin kw + ig sin kw 

which is real only when sin kw is zero. This implies that cos kw is ±1. 
Thus if | fci | < 1, G (jw) cannot encircle the +1 point. (This result is 
in agreement with Assertion 2 of Ref. 5, to which this problem corre- 
sponds if gr = 0. It is intuitive that adding losses to a lossless system will 
enhance stability.) 

To show the necessity of the open-loop stable requirement, note that 
g can be either positive or negative. From the Pontryagin Criterion, we 
see that (27) is then stable or unstable, respectively, and that for 
| kt | < 1 the closed-loop system is stable or unstable, respectively. 
However, in either case there are no encirclements of the critical point 
by the open-loop gain. 
Example 3: Let 
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2 a « s ' 

G(s) = -^ b n ^ 

(28) 
p 

z c y(s) = ^ — Cp ^o. 
E<te' 

In this final example we look at how the a, b, c, and d coefficients of (28) 
enter into the R t (s) polynomials of (6) and into the Q„(e~ z ) polynomials 
of Theorem 1 and its corollaries. We show how the Assumptions A2 
and A3 affect whether the system satisfies Corollaries 1 and 2, and 
thereby provide examples of such systems. Using (3) and (28), (5) 
becomes 

F(s) = s n+p \b n (d p + c„) - 2a n d p e- ks + b n {d v - cje -8 *'} 

+ s-^-'JIXdp-! + c p _,)b n + b n ^(d p + c p ) 

— 2[a n . l d p + Mp-Je - ** 

+ [&„(</„_, - c,,.,) + &„_,(<*, - c p )]e- 2k '\ 

+ "Es"I E W* + c) - 26-*' 2 M< 

m-0 i+i=m j + i-m 

+ e" 2 " E W*-cO). (29) 

I + i = m 

First we investigate the zeros of the coefficient of s n+p in (29). [This 
coefficient corresponds to Q„(e~ : ) in Corollary 2.) For simplicity let 
e" = z. This maps the left-half s-plane into the unit circle in the z plane. 

If (/„ = c p , then the coefficient of s n+p has zeros whenever b n c p = 
a„Cp2~ l . If a n were zero, the coefficient in question would become constant, 
which satisfies the conditions of Corollary 2. If a n is not zero, then the 
condition under discussion simplifies to 2" 1 = b n /a„ . All solutions of 
this will satisfy | z \ < 1 if | a„ \ < \ b„ \ . Hence all zeros of the coefficient 
of s n+ " in (29) will lie in the left-half plane if | a n | < \b n \ . This is 
intuitively appealing since this requires that G(s) have less than unity 
gain at large frequencies (as required by Assumption A2). 

On the other hand, if d p 9^ c p , then the zeros of interest are solutions of 

<*-'>'- W^' + ^ = °- (30) 
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It is well-known 6 that solutions of (30) (for z _1 ) have magnitude less 
than unity if and only if the following three conditions are met. 

i+-& < 1 (31a) 

dp — c p 

d v + c p _ 2d p 2aJ p ,„.,. 

1 + d v - c p - d p - c p > b n (d p - c p ) {61b) 

24. > _ 2oA (31c) 



rf P — c p b„(d p — c P ) 

Thus conditions (31) are NAS for | z \ > 1. Condition (31a) requires that 
d p and c p have opposite sign. Since this corresponds to terminating the 
line in a negative conductance at high frequencies [i.e., lim,_ >00 Y(s) < 0], 
we reject this case. If both d p and c p are nonzero, and have the same 
sign, (31a) is violated. If d p ^ 0, (31b) and (31c) together require 
| a n | < | 6„ | as before. The remaining possibility is that d p = 0. This 
is a reasonable physical assumption; in fact, it is required by Assump- 
tion A3. Invoking Assumptions A2 and A3, the coefficient in question 
now becomes 6„c p (l — e~ 2k ") which has an infinity of purely imaginary 
zeros, and this example no longer satisfies Corollary 2. 
To see if it satisfies Corollary 1, rewrite (29) as 

F(s) = s n+p -M(s&nC p + 6 n -,c p + 6„c p _,)(l - e" 8 **) - 2a„d p -ie - *' 
+ b n d p ^(l +e- 2k ')\ 

+ " + E s m { E IhWi + c.) - 2e-**a,di + e^'b^d, - c,)]} . (32) 



We complete this example by finding conditions under which the coeffi- 
cient of s n+p_I in (32) satisfies the conditions of Corollary 1. This coeffi- 
cient can be written as 



ks ^ + b n . x c p + 6 n c p _, + Mp-, - 2a„d p _ 1 e "' 



- e- 2kl 



(ks^ + b„-ic p + &„c p - t - bj^ij- 
Let w = ks. Then (33) becomes 

A; L K c p c p J 

_ _„ 6„c p f . fcb n _! fec,_i fed,-, "! a. 



(33) 



(34) 
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Zeros of (34) are given by the solutions of (35) 

(10 + a) + 77" - (w + &)e~ 2w = (35) 

where 



'-■ft? 



c p c v 

_i_ c p-i »p-i 



7 = -2k a f^. 
b n c p 

We assume that a > /3 since a — = 2d p _i/c p fc which is positive by 
Assumption A3. 

It is also reasonable to assume a > 0, since b„- l /b„ must exceed zero 
for the denominator of G'( s ) to be strictly Hurwitz,* and since c p -i/c p 
less than zero would imply zeros of Y(s) in the right-half plane. These 
considerations also imply that \ a \ > | /3 | . Using these assumptions 
(i.e., a > ft a > 0, | a | > | |) it follows that 

| to + a | > | IP + /3 | 

for Re (w) ^ 0. Evaluating the magnitude of (35) on Re (w) ^ yields 

| w -\- a + ye'" — (iv + 0)e~ 2w \ 

^ | w -f a | — | 7 | | e -1 " | — | 10 + /3 I I e~ 2w | 

> | 10 + a | — I 7 I — |i0 + /3|. 

If lim.-.ao G(s) = 0, then a„ = and 7 = 0. This means that (35), and 
hence (33) and (34), have no zeros in the right-half plane and thus 
Corollary 1 is applicable to this problem. 

VI. CONCLUSIONS 

It has been shown that the time-honored technique of determining 
the existence of unstable poles of a closed-loop gain by counting encircle- 
ments of the critical point of the open-loop gain along a finite segment 
of the imaginary axis remains valid for a large class of distributed 
parameter systems of practical importance for which the open-loop gain 
approaches zero for large frequencies. Existing limitations of the Michai- 



* A well-known necessary condition for a polynomial to be Hurwitz is that all 
coefficients have the same sign (see, for example, Ref. 6, p. 281). 
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lov Criterion have been removed so as to include physical systems of 
lossless transmission lines. 
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APPENDIX 

a.i Pontryagin's Criterion 

Let h(z, t) be a polynomial with complex coefficients in the two 
complex variables z and t. Pontryagin 1 has developed necessary and 
sufficient conditions that the function H{z) = h(z, e*) have zeros with 
only negative real parts. We now present one of Pontryagin's main 
results. Let r and s be the degrees of the polynomial h(z, t) with respect 
to z and t. Then the principal term of h(z, t) is the term containing the 
product z T V. Pontryagin showed that if h(z, t) does not contain the 
principal term, then H (z) has an infinity of zeros with arbitrarily large 
positive real parts. 

Let p(-) and q(-) be real-valued functions of a real variable. We say 
that the zeros of these two functions alternate if: {i) they have no 
common zeros, (u) they have only simple zeros, and (in) between every 
two zeros of one of these functions there exists at least one zero of the 
other. The result of Ref. 1 which will be used in the present study is: 

Let h(z, t) be a polynomial with the principal term and H(iy) = 
F(y) + iG(y) where F(y) and G(y) take on real values whenever y is 
real. If all zeros of the function H(z) have negative real parts, then all 
zeros of F(ij) and G(y) are real, alternate, and 

0'(y)F(y) - F'(y)G{y) > 

where superscript prime denotes the derivative. In order that all zeros 
of H (z) have negative real parts, it is sufficient that all zeros of F(y) 
and G(y) are real and alternate and that G'(ij)F(y) - F'{y)G{y) be 
positive for some y. 

a.2 Proof of Theorem 1 

Choose a real number R' > 1 such that Re (z) > and F k (z) = 
implies \ z \ < R'. Define the set 6, 

d = {z|Re(z) >0 and \z\> R'.}. 
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Let 

D k = inf | F k (z) | > 

tie 
and 

M k = sup 2| a ti |. 

i-0.---.k-l t = l 

If M k = 0, then F(z) = F k (z) and the theorem is trivially true. If 
M fc > 0, then for all z t d the following inequalities hold. 

k i 



F(z) | ^ | ***•,(«) | - £ | *'Q,(0 | 

1-0 

^ I ** I J>* - E I *' I E I «■-.• 

1-0 t=l 

1-0 t=l 

^ £* | 2 |* - M k £ J 2'' I 

£ Z>* I z I* - M t 



f-0 

lit 1 

2 | — 1 

|2 I - 1 



2 |* [D k (\ z 1 - 1) - M k ] + M t 
1*1-1 

[D» I • I - (Di + MO] + tti 



- |« I - 1 

which is positive for | 2 | > 1 + M k /D k = R k . Hence the magnitude of 
all zeros of F(z) must be bounded by R, the larger of the numbers R' 
and R k . The theorem follows at once by noting that F(z) is analytic 
and that an analytic function has at most a finite number of zeros in 
any finite region. 

A.3 Proof of Theorem 2 

We now wish to derive an expression for the number of rhp zeros 
of F(z). We choose a contour r varying along the imaginary axis from 
— y to y (call this portion o>) where y ^ R k of Theorem 1 and close it 
by a contour C outside the semicircle of radius R of Theorem 1. Let 
F(z) and \p(z) be as defined in (15) and (16), respectively. We choose 
contour C (and increase y, if necessary) so that | yf/{z) \ < 1 along C. 
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Let N be the number of zeros of F(z) inside T and let 

A r (F(js)) 

be the net change in arg (F(z)) along T. Then N, the number of zeros 
of F(z) enclosed by T (assuming counterclockwise travel), is given by 



N = ^ A r (F(z)) = f A u (F(z)) + ^ A c (F(a)) 



= ;p A„(F(z)) + A c (/) + X(F k (z)) + A c (l + *(z)) 

since F(z) has no zeros on T. Since 1 + <A(z) does not wind around the 
origin (its real part being always positive), A c (l + ^(z)) = 0. Since 
A c (z k ) = fc/2 we have thus proven Theorem 2. 

A.4 Sufficient conditions that Q n does not wind around the origin 
Let 



= mm r, 

j-2. •••.m 



Then 



e — Re [Q n ] = e rkz + £ ~ <T <r '" rl)l -cos r t y 

a ln ,=2 &ln 



Re [QJ I ^ <*'" - £ 



Ofa 



I «m 

Hence | Re [Q n ] | > f or x > \/r k In a where 



«« E 



a ln 



Thus we need only consider a rectangle defined by < x ^ In a/r k , 
£ y < P. [Note that any contour C will work if a < 1 which is the 
case if assumption (11) is used.] Q n (e~°) does not wind around the origin 
for any contour with x > In a/r k , since Re [Q„] does not change sign. 
If E7-i o7»" 5^ for < a; < In a/r* then Re [Q n ] does not change sign 
for y any integer multiple of P, and < x < In a/r* . A suitable contour 
could then consist of a horizontal line at y = KP from x = to x = 
In a/r k , where K is an integer large enough so that KP > R of Theorem 1. 
The rest of the contour in the first quadrant could be semicircular. The 
contour is completed in the fourth quadrant by the mirror image of the 
first quadrant. Since Re [Q n ] 9* along this contour, Q n [e~ z ] does not 
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wind around the origin. (Note that a in > 0, j = 1, — , m is sufficient 
to satisfy the conditions of this paragraph.) Furthermore, since Re [Q n ] 
is even in y, A c (Q n (e~ ! )) = along the contour chosen. 

This result can be extended to include the case where Re [Q n ] has 
simple zeros on y = KP. In this case, use semicircular indentations 
around such points, in the direction to have Im [Q„] > 0, in both first 
and fourth quadrants. (Hence the contour ceases to be symmetrical 
about the real axis.) Thus, along the deformed horizontal lines, the 
graph of Q„(e~') remains in the upper-half plane. Along the semicircular 
portion it remains in either the right- or left-half plane. Hence no 
encirclements of the origin are possible and again A c [Q„(e~*)] = 0. 
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